LAGUERRE AND MEIXNER SYMMETRIC FUNCTIONS, 
AND INFINITE-DIMENSIONAL DIFFUSION PROCESSES 



GRIGORI OLSHANSKI 

Abstract. The Laguerre symmetric functions introduced in the note are indexed 
by arbitrary partitions and depend on two continuous parameters. The top degree 
homogeneous component of every Laguerre symmetric function coincides with the 
Schur function with the same index. Thus, the Laguerre symmetric functions 
form a two-parameter family of inhomogeneous bases in the algebra of symmet- 
ric functions. These new symmetric functions are obtained from the A-variate 
symmetric polynomials of the same name by a procedure of analytic continuation. 
The Laguerre symmetric functions are eigenvectors of a second order differential 
operator, which depends on the same two parameters and serves as the infinitesi- 
mal generator of an infinite-dimensional diffusion process X(t). The process X(t) 
admits approximation by some jump processes related to one more new family of 
symmetric functions, the Meixner symmetric functions. 

In equilibrium, the process X(t) can be interpreted as a time-dependent point 
process on the punctured real line R \ {0}, and the point configurations may be 
interpreted as doubly infinite collections of particles of two opposite charges and 
log-gas-type interaction. The dynamical correlation functions of the equilibrium 
process have determinantal form: they are given by minors of the so-called ex- 
tended Whittaker kernel, introduced earlier in a paper by Borodin and the author. 
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1. Introduction 

1.1. Preface. The present note is a research announcement; the detailed exposi- 
tion will appear elsewhere. The goal of the work is twofold: (1) to introduce new 
bases {£„} and {OJt^} in the algebra A of symmetric functions, and (2) to construct 
a diffusion process X(t) in an infinite-dimensional cone Q. The two subjects are 
interrelated: The algebra A serves as the algebra of "polynomial observables" on 
the cone Q, and the basis elements £„ e A are the eigenfunctions of the infinites- 
imal generator of the process X(t). As for the basis {971^}, its elements are the 
eigenfunctions of the infinitesimal generator of an auxiliary jump process. 

The basis elements 2, v G A are called the Laguerre symmetric functions. They are 
indexed by arbitrary partitions v and depend on two continuous parameters. The 
diffusion process X(t) depends on the same parameters. It possesses an invariant 
symmetrizing measure W, which also serves as the orthogonality measure for the 
Laguerre symmetric functions. 

1.2. Finite-dimensional counterparts. All the basic objects - - the Laguerre 
symmetric functions, the cone Q, the probability distribution W on Q, and the 
diffusion X(t) on Q — have finite-dimensional counterparts; I will describe them 
briefly. 

• In dimension 1, these are the classical Laguerre orthogonal polynomials on the 
half-line x > with the weight measure x b ~ 1 e~ x dx (here b > is a parameter), and 
the diffusion is generated by the associated ordinary differential operator 

d 2 . d 

the Laguerre polynomials are its eigenfunctions. 

• In dimension N = 2, 3, . . . , we deal with the algebra Ajv of symmetric polyno- 
mials in N variables xi, . . . ,xn- Such polynomials are viewed as functions on the 
A-dimensional cone 



d N = {x:= (x 1 , . . . , x N ) : x\ > ■ ■ ■ > x N > 0} C 



aJV 
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A relevant basis in is formed by the A-variate symmetric Laguerre polynomials 
[Ma87j . |La91cj . which are orthogonal with respect to the measure on Q N with the 
density 

(x 1 ...x N ) b ~ 1 e-^-- XN - H (x t -x 3 f. 

l<i<j<N 

Assuming x to be in the interior of the cone, one can interpret x as a collection of N 
indistinguishable particles on the half-line M>o; then the above measure determines 
an ensemble of random particle configurations, called the A-particle Laguerre en- 
semble. Again, there exists an associated diffusion process Xpr(t) with state space 
Qn. In the interior of the cone, X^(t) may be interpreted as a random evolution of 
N interacting particles with pairwise repulsion. One may call X^(t) the A-particle 
dynamical Laguerre ensemble. 

• Note also that there exists a lattice version of the Laguerre ensemble, related to 
natural discrete analogs of the Laguerre polynomials — the Meixner polynomials. 
There exists also an associated Markov dynamics, which is a Markov jump process. 
In the simplest case N — 1 this is a birth-death process with linear jump rates. 

1.3. Analytic continuation. In the literature, there exist many models of such a 
kind, continuous and discrete, associated with various systems of orthogonal poly- 
nomials. A general recipe for building infinite-dimensional analogs of such models, 
often applied in Random Matrix Theory, is to use a large- A limit transition (see, 
e.g., the survey paper |KT10cj ). However, in the present work, a different approach 
is applied. In short, its main idea can be formulated as extrapolation into com- 
plex domain via analytic continuation with respect to two parameters, the number of 
particles N = 1, 2, 3, . . . and the additional parameter b > 0. Surprisingly enough, 
although parameters N and b are of a very different origin, they can be treated on 
equal grounds. 

In comparison with the existing approaches to random-matrix-type dynamical 
models in infinite dimension (see [SpBT] , [Os09] . |KT09j . jKTlOaj . |KT10b] K our 
approach is to great extent more algebraic. 

1.4. Point processes. In the equilibrium state corresponding to the stationary 
distribution W, the process X(t) can be interpreted as a time-dependent point 
process X stat (t). This interpretation relies on the fact that the invariant measure 
W is supported by a dense Borel subset Q' C Q admitting a natural realization as 
a space of infinite particle configuration on the punctured real line R* := R \ {0}. 
It turns out that the time- dependent point process X stat (t) is determinantal: its 
dynamical correlation functions are given by minors of a kernel K(s, u; t, v) on (M x 
M*) x (1 x 1*) — the so-called extended Whittaker kernel, which initially appeared 
in [BO06aj . 

Note that the A-particle dynamical ensembles Xjy(i) live on the half-line R>o 
while the particle configurations of the process X stat (t) live on the punctured real 
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line R*, which is the union of two half-lines. Because of this duplication effect, the 
claim that X(t) cannot be related to the processes X^(t) through a large- N limit 
transition becomes intuitively evident. 

1.5. Lattice approximation. Although X(t) does not arise from a large- N limit, 
it admits a lattice approximation. Namely, X(t) can be obtained as a scaling limit 
of some jump processes depending on an additional parameter £ G (0, 1), as £ X 1. 
These jump processes were studied in detail in [BO06aJ. Their state space is the 
set of all partitions. This countable set also can be realized as a set of particle 
configurations on the lattice Z + | of half-integers; here the number of particles is 
finite but not restricted. The second basis {97^} C A mentioned above just arises 
in connection with these jump processes. The elements VJl u are called the Meixner 
symmetric functions; they depend on three parameters: to the two parameters of 
the Laguerre symmetric functions one adds the third parameter £ G (0, 1). In the 
limit as £ y 1, one has D3l u — > £„, similarly to the approximation of the classical 
Laguerre polynomials by the Meixner polynomials. 

1.6. Concluding remarks. The results announced in the present note continue 
those of jBUOO] . |BU06a] . |BO06bj . |BU09] . and all these works have a connection 
with the asymptotic representation theory of the symmetric groups. It is interesting 
to compare the results of these papers and the present note with those of the pa- 
pers |BO05aj . [BO05bJ, |BO10j . which are related to representations of the unitary 
groups. Although the construction of a Markov dynamics in [BOlOj relies on a dif- 
ferent approach, |BO10t Appendix] also exploits the idea of analytic continuation. 
Note that in the context of the unitary groups, the Laguerre polynomials are (to 
some extent) replaced by the Jacobi polynomials. 

Finally, I would like to note that the representation theory of reductive groups and 
Lie algebras also affords examples in which finite-dimensional objects ciriSG clS db de- 
generation of infinite-dimensional ones, and, conversely, infinite-dimensional objects 
may be reconstructed from finite-dimensional ones through analytic continuation 
in parameters. For instance, the principal series representations or highest weight 
modules may be viewed as analytic continuation of the irreducible finite-dimensional 
representations. 

2. The one-particle case 

There is a well-known relationship between systems {<p n ; n = 0, 1, 2, . . . } of orthog- 
onal polynomials of hypergeometric type on R and some one-dimensional Markov 
processes X{t) (see, e.g., |Sch00] ). Namely, the state space of X{t) is a closed subset 
I C R — the support of the orthogonality measure of {<p n }, and the infinitesimal 
generator of X(t) is a second order differential (or difference) operator D, such that 



D(p n = -/i„0„, [i n > 0. 
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The orthogonality measure w for the polynomials <p n serves as an invariant and 
symmetrizing measure of X(t). The transition function of X(t) has the form 

The simplest examples are provided by the classical Jacobi, Laguerre, and Hermite 
polynomials: 

• Jacobi polynomials: I is the closed interval [—1, 1], w has density (1 —x) a ^ 1 (l + 
x) b ~ x with parameters a, b > 0, and 

D=(l-x 2 )^ + [b-a-(a + b)x]± 

• Laguerre polynomials: I is the closed half-line [0,+oo), w has density x e~ x 
with parameter b > 0, and 

^ 2 /, \ d 
D = x— - + [b- x)—. 

dx z ax 

• Hermite polynomials: I is the whole real line, w has density e~ x2 ^ 2 , and 

D _ d 2 

dx 2 dx 

In the Hermite case, X(t) is the Ornstein-Uhlenbeck process. In the Laguerre 
case, X{t) is closely related to a squared Bessel process (see, e.g. |Eie83j ). 

3. The TV-particle case 

Here we recall a well-known construction providing a multidimensional general- 
ization of the above picture. 

Fix N = 1,2,3,... . Instead of univariate polynomials we will deal with sym- 
metric polynomials in N variables x\, . . . ,x^. Denote by A^r the algebra of such 
polynomials (the base field is R or C depending on convenience). The interval / is 
replaced by the subset 

J OTd : = {(^l) • • • 3 x n) C I n : Xi > ■ ■ ■ > x N }, 

and we regard A N as an algebra of functions on I^ d . Let v = (z/ 1; . . . , v N ) range 
over the set of partitions of length t{y) < N. We set 

<I> v \n{xx,...,x n ) = —^ 3.1) 

V N [xi,...,x N ) 

where the determinant is of order N and Vn is the Vandermonde, 

V N = V N (x u ...,x N )= J [ (xi - Xj). 

l<i<j<N 
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The <p u \pf's are symmetric polynomials forming a basis in Ajy. Moreover, it is readily 
verified that they are pairwise orthogonal with respect to the measure 



N 



Wn{(1xi . . . dxx) = (Vn) 2 Y\ w(da 



(3.2) 



on j£ d . 



This construction seems to be well known; see e.g. Lassalle's papers [La91a] , 
[La91bj . |La91cj . Formula (13. ip is similar to the classical expression for the Schur 
symmetric polynomials (the Schur polynomials appear if one substitutes 4> n {x) = x n \ 
they are not orthogonal polynomials though). 

Further, the analog of D is the operator 



D N := V N l o (D Xl + ■■■ + D XN ) o V N + d N l 
where D Xi stands for a copy of D acting on variable Xi and 

u>n = fio + " • + Hn-i- 
If D is a differential operator, 

D = A {X )^ + B {X )± 
then D]\r is a partial differential operator, 



(3.3) 



N 



D 



N 



A ^ 



i=l 



dX 2 : 



+ 



B(xi) + 



2A(xi 



X j 




Due to the special choice of dpj, the constant term in vanishes. 

Although the coefficients in front of the first derivatives have singularities along 



the hyperplanes Xi 



x 



the operator D]y is applicable to symmetric polynomials 



and preserves the space A^. The polynomials 4> u \n are eigenfunctions of D 



N- 



N 



D 



N<Pu\N 



>v\N- 



J=l 



Finally, D N serves as the (pre-)generator of a Markov process X N (t) on J^ d with 
invariant symmetrizing measure w^. 

The first example of such a process X N {t) has been investigated in |Dy62| ; it 
corresponds to the system of Hermite polynomials. As shown in that paper, X^it) is 
obtained from a matrix-valued Ornstein-Uhlenbeck process through the projection 
onto the matrix eigenvalues. One may say that X^if) is the radial part of this 
matrix-valued Markov process. 

In the present work we focus on the Laguerre case. 



LAGUERRE AND MEIXNER SYMMETRIC FUNCTIONS 



7 



4. The Laguerre symmetric functions 
Let ei, ei, ■ ■ ■ denote the elementary symmetric polynomials, 

C\ } j ^ii ^2 ^ ^ XjXj , 63 ^ ^ X{XjX k , 
i i<j i<j<k 

and so on. Here it is tacitly assumed that the indices range over {1, . . . , N}, where 
N is the number of variables. As well known, the algebra Ajv of iV-variate symmetric 
polynomials is isomorphic to the algebra of ordinary polynomials in e\, . . . , e^. 

Our first step is to make a change of variables: take {e 1 , . . . , e^} as new (formal) 
variables instead of natural coordinates x\, . . . , xn- 

Theorem 4.1. The N-variate Laguerre operator D N : A N — > A N can be rewritten 
as the following second order partial differential operator in variables e\, . . . , e^: 

N f>2 N o 

D N = 2^ A mn +2^ B n 



where 



and 



de m de n ^ " „ ' 



min(ra,n) — 1 

A m n = ^ (m + n-1 - 2k)e m+n - 1 - k e k 

k=0 



B n = -ne n + (N - n + l)(iV + b - n)e„_i 
to^/i t/ie convention that eo = 1 and = /or k > N. 

The next step is to replace A N by the algebra A of symmetric functions. For 
our purpose, it is convenient to define A as the algebra of polynomials in countably 
many formal commuting variables e±,e2, ■ ■ ■ , which are assumed to be algebraically 
independent. 

For N = 1,2, ... , let Jn C A denote the ideal generated by elements e k with 
k > N. The quotient algebra A/Jn is naturally isomorphic to A^, so we get a 
canonical algebra homomorphism hn '■ A — > An, which we call the Nth. truncation 
map. 

Definition 4.2. Let z and z' be complex parameters. Consider the formal second 
order differential operator in countably many variables e\, ■ ■ ■ , obtained from the 
A^-variate Laguerre operator Dn by removing the relations eN+i = eN+2 = • • • = 0, 
dropping the restriction m,n < N, and replacing the factor (N — n + 1)(AT + 6 — n) 
in the definition of coefficient B n by [z — n + l)(z' — n + 1): 



m,n=l d6md6n ™ 9e «' 



8 



GRIGORI OLSHANSKI 



where A mn is given by exactly the same formula as above, 




Observe that D is correctly defined as an operator A — > A; we call it the Laguerre 
operator in A. 

If z = N and z' = N+b—1, then D preserves the ideal Jn C A and hence factorizes 
to an operator in the quotient A/ J N = A N ; the resulting operator coincides with D N . 
Note that this property, combined with the polynomial dependence on parameters 
z, z', determines the operator uniquely. In this sense, D may be viewed as the result 
of analytic continuation (or extrapolation) of the A^-variate Laguerre operators 
with respect to parameters A" and b. 

Let {L n } denote the system of monic Laguerre polynomials with parameter b > 
(see, e.g., |KS96j ). Recall that, in our notation, the weight function is x b ~ 1 e~ x . 
Next, let {L u \ Nb } denote the system of A^-variate symmetric Laguerre polynomials 
defined in accordance with the determinantal formula (13. ip . We are going to define 
elements of A that may be viewed as are analogs of the polynomials L u \ Njb . To do 
this we apply the same principle of analytic continuation in A" and b as we have 
employed in the definition of £>. 

Theorem 4.3. For an arbitrary partition v = (u±, V2, ■ ■ ■ ), there exists a unique 
element £„ e A[z, z'\ :— A <g> C[z, z'\, such that for any natural number N > iiv) 
and any b > 0, 



Here ti^ : A — > An is the A^th truncation map introduced above. 

Definition 4.4. We call the elements £ u the Laguerre symmetric functions. 

Recall the definition of the Schur symmetric functions: these are elements of A 
indexed by arbitrary partitions v and expressed through the generators e n by the 
following formula (the Nagelsbach-Kostka formula, see pvla95j) 



here v ' stands for the partition given by transposing the Young diagram correspond- 
ing to z/, and the order of the determinant is an arbitrary integer > i{y'\ As well 
known, the Schur functions form a distinguished homogeneous basis in A, and 



Let us explain the notation used in the next theorem. We identify partitions and 
the corresponding Young diagrams. Given a couple \i C v of Young diagrams, we 




S u — det[e„/_j+j]; 
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denote by dimzz/yU the number of standard Young tableaux of the skew shape vf/i. 
The symbol □ e u/jj, denotes a box in v/p,, and c(D) denotes its content, equal to 
the difference j — % of the column number j and the row number % of the box. 

Theorem 4.5. TTie expansion of the Laguerre symmetric function £ u in the basis 
of the Schur symmetric functions has the form 

where 

l7(,,^,,/) = (-i)Hhh^W^ jj {z + ci nw + c(n))). 

U I \t i \)- ngl/ / M 

Since C{y, v; z, z') — 1, the top homogeneous component of £„ is equal to S u : 

£v — S v + lower degree terms. 

It follows that the Laguerre symmetric functions with any fixed values of parameters 
z and z' form a basis in A. 

For the empty diagram corresponding to the zero partition, v = 0, we have 
£j = S = 1. This is the only case when Z v and S v coincide: for v ^ 0, SL V is an 
inhomogeneous element, so that the basis {£ u } of Laguerre symmetric functions is 
an example of inhomogeneous basis. In this respect, it differs from other bases in A 
that are commonly used in algebraic combinatorics. 

A box □ in a Young diagram v is said to be a corner box if the shape v \ □ 
obtained by removing □ from v is again a Young diagram. Let v~ denote the set 
of all corner boxes in v. For instance, if v = (3, 2, 2) then v~ comprises two corner 
boxes, □ = (1, 3) and □ = (3, 2). 

Theorem 4.6. The action of D on the Schur functions is given by 

®S V = -\u\S v +J2( Z + c (°))( / + c ( Q ))^\n 
nev- 

Theorem 4.7. The Laguerre symmetric functions are eigenvectors of the Laguerre 
operator D with the same values of parameters (z,z f ). More precisely, 

D£ v = -\v\2 v . 

5. Formal orthogonality 

Definition 5.1. For any fixed (z, z') e C 2 , introduce the formal moment functional 
tp : A — > C by setting 

^(1) = 1, ^(3)/)=0 for any fe A, 

where the Laguerre operator D is taken with the same values of the parameters as 

1>- 
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The definition is correct by virtue of the last theorem. Indeed, it implies that 
that the range of D is the span of the Laguerre functions £ u with v ^ 0, while the 
vector 1 = £ e A is transversal to this span. Note that ip depends polynomially 
on the parameters z, z' . 

Theorem 5.2. For any Young diagram v , 

^(S u ) = H(z + c(D))(z' + c(D))- 

This formula provides an alternative (equivalent) way of introducing the moment 
functional. 

Theorem 5.3. For any two Young diagrams v and \x, 

V = S vfl ]J(z + c{U)){z' + c(D)), (5.1) 

where 5 Vfl is Kronecker's delta. 

This result shows that the Laguerre symmetric functions are pairwise orthogonal 
with respect to the inner product in the space A defined by 

(f,g) = {Lg)z, Z '--=ip(f9), f,geA. (5.2) 

Obviously, the inner product is nondegenerate if (and only if) z and z' are not 
integers, for then the product in the right-hand side of (15. ip never vanishes. 

Remark 5.4. Assume z and z' are not integers. Then the Laguerre symmetric 
function £ u is characterized by the following two properties: 

(1) £, u differs from the Schur symmetric function S u by lower degree terms; 

(2) Z v is orthogonal, with respect to inner product (15. 2p . to all elements of A of 
lower degree (that is, of degree strictly less than deg£y = \v\). 

Alternatively, without any assumption on the parameters, 2, u is characterized by 
(1) together with the following property replacing (2): 
(2') 2j v is an eigenfunction of the Laguerre operator D. 

6. The orthogonality measure for the Laguerre symmetric 

functions 

Definition 6.1. By the Thoma cone we mean the subset f2 C x E2° x M + 
consisting of triples u> = (a, f3, 5), where 

a = {a 1 >a 2 > ■■■>0), (3 = (ft > > • • • > 0), 5>0 

and 

oo 

1=1 



/ dim i 
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The Thoma simplex Vt is the subset of fl determined by the additional condition 
5 = 1. 

Both & and f2 are closed subsets of the product space M 2oo+1 := M°° x M°° x M 
equipped with the product topology; Q is compact and fl is locally compact; fl is 
precisely the cone with base fl. 

Let Pun(Q) denote the space of continuous functions on the Thoma cone; this is 
an algebra under pointwise multiplication. We embed A into Fun(f2) by setting 

oo oo 

l + £*rf"->e*n£fj. 

n=l i=l 

where t is an auxiliary formal variable and 

oo 
i=l 

More precisely, the correspondence, which is written in terms of the generating series 
for {e„}, turns each e n into a function e n (u) on the Thoma cone. This function is 
continuous. The correspondence e n i— >■ e n ( • ) is extended to the whole algebra A by 
multiplicativity. In this way we get an algebra morphism A — > Pun(fi), which is an 
embedding. For any / G A we denote the corresponding continuous function on fl 

by/H. 

Equivalently, in terms of another system of generators of the algebra A, the New- 
ton power sums p n , the embedding A — > Fun(f2) can be defined by setting 

oo oo 

Pl (u)=5, Pn (o;) = ^<-^(-Ar forn>2. 

i=\ i=l 

In the algebra A, there is a distinguished involutive automorphism, which is de- 
fined on the generators p n as p n h- >■ (— l)™^ 1 ^. Under this automorphism, S u h- >■ S^/. 
The above formula shows that in the realization A C Fun(Q), this automorphism 
amounts to transposition a -H- /3. 

Definition 6.2. Let us say that (z, z') G C 2 is admissible if both z and z' are 
nonzero and 

Y[(z + c(D))(z' + c(D))>0 

for any Young diagram v. 

The family of admissible couples (z, z') splits into the union of the three subfam- 
ilies: 

• The principal series: both z and z' are no meal, z' = z. 

• The complementary series: both z and z' are real and are contained inside an 
open interval (m, m + 1) with m G Z. 



12 



GRIGORI OLSHANSKI 



• The degenerate series: (z, z') = ±(iV, A" + b — 1) or (V, z) = ±(iV, AT + ft — 1), 
where iV = 1, 2, . . . and b > 0. 

Theorem 6.3. Assume (z, z') is admissible. There exists a probability measure 
W = W z ,z' on the Thoma cone Q such that all functions from A are W -integrable and 
the formal moment functional ip with parameters (z, z') coincides with expectation 
with respect to W : 

^(/)= ff(uj)W(du) for all fe A. 
Jo, 

These properties determine W uniquely. Moreover, the functions from A are dense 
in the Hilbert space L 2 (Q, W). 

By virtue of Theorem 15.31 this implies that the measure W with admissible pa- 
rameters (z, z') is the orthogonality measure for the Laguerre symmetric functions. 

In the case when z — N — 1, 2, . . . and z' = N + b — 1 with b > we recover the 
Laguerre measures wn on the A^-dimensional cone 

'■= {(xi > ■ ■ ■ xn > 0)}. 

Here we embed Qn into the Thoma cone by setting «j = Xi for % = 1, . . . ,N and 
5 = ^Xi (so that all remaining a- and /3-coordinates are equal to zero). 
An immediate consequence of the theorem is 

Corollary 6.4. Assume that the couple (z,z') is admissible and nondegenerate, i.e. 
belongs to the principal or complementary series. Then the Laguerre symmetric 
functions, viewed as functions on the Thoma cone, form an orthogonal basis in the 
Hilbert space L 2 (Q, W z z >). 

7. Properties of the orthogonality measure 

Here are some properties of the measures W ZtZ / with admissible parameters (z, z')\ 

• W Z)Z i are obtained from the so-called z-measures on the Thoma simplex by a 
simple integral transform along the rays of the Thoma cone. See |BO00t Section 5]. 

• W z ,z> does not change under transposition z -H- z' . Within this symmetry 
relation, the measures corresponding to different couples of parameters are pairwise 
disjoint: this follows from |KO V04] . 

• The involutive map a -H- f3 of the Thoma cone transforms W z z > to W- z - z >. 

• W e z> is supported by the subset of the Thoma cone formed by those triples 
lo = (a, (3, S) for which 5 exactly equals J2( a i + A)> i- e - 7 = 0. This subset is Borel 
and everywhere dense in Q. (Note that 7 is not a continuous function in w, it is 
only lower semicontinuous.) 

Theorem 7.1. If (z,z') is in the principal or complementary series, then the topo- 
logical support ofW z ^ z i is the whole space Q. 
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(The topological support of a measure is the smallest closed subset supporting 
the measure.) 

The theorem implies that a nonzero continuous function on Q, which is square 
integrable with respect to W z , z /, is a nonzero vector in the Hilbert space L 2 (fl, W Z)Z /). 
In particular, the natural map assigning to a bounded continuous function a vector 
in this Hilbert space is an embedding. This assertion fails in the case when (z, z') 
belongs to the degenerate series. 

8. The Laguerre diffusion process on the Thoma cone 

In this section (z, z') is a fixed couple of parameters from the principal or comple- 
mentary series. 

Recall that we may regard A as a dense subspace in the Hilbert space H : = 
L 2 (Q, W z>z >) and then {£„} becomes an orthogonal basis in H. The Laguerre oper- 
ator Q : A — > A is diagonalized in this basis, 

= -\v\2», 

so that the eigenvalues of D are 0, —1, —2, . . . , where has multiplicity 1 and corre- 
sponds to the basis vector £ = 1. It follows that D generates a strongly continuous 
semigroup P(t) of contractive selfadjoint operators in H, 

P{t)Z u = e -1 " 1 *^, t > 0. 

Theorem 8.1. The semigroup P(t) is a conservative Markov I? -semigroup 

By definition, this means that P{t)l = 1 and P(t) preserves the cone of nonnega- 
tive functions in H. The first claim is obvious, the second one is nontrivial; its proof 
relies on the approximation by some jump Markov processes, see Section HH 

The next claim says that P(t) is actually a Feller semigroup (in one of the versions 
of this property). Let C(ft) C Fun(fi) be the space of bounded continuous functions 
on Q with the supremum norm and Cq(Q) C C(fl) be its closed subspace formed 
by the functions vanishing at infinity. Both C(Q) and Cq(Q) are Banach spaces 
contained in H, but Cq(Q) is separable while C(Q) is not. 

Theorem 8.2. The semigroup P(t) preserves Cq(Q) and induces a strongly contin- 
uous contractive semigroup in this Banach space. 

One of the ingredients of the proof is separation of variables described in the next 
section. 

The theorem implies that P(t) gives rise to a Markov process X(t) on Q with 
cadlag sample trajectories. Actually, more can be said: 

Theorem 8.3. With probability one, the trajectories of X(t) are continuous. 

By the very construction of the Markov process X(t), the probability measure 
W z>z i is its invariant and symmetrizing measure. 
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9. Separation of variables 

Extend the algebra A by allowing division by e± (that is, localize over the multi- 
plicative semigroup generated by ei) and denote the resulting algebra by A ext : 

A ext :=C[e 1 ,er 1 ;e 2 ,e 3) ...]3A. 

Since the coefficients of the Laguerre differential operator D : A — > A (Definition 14.21 
are polynomials in variables e n , D can be extended to an operator D ext : A ext — > A cxt . 
Set " 

r := ei, e° := e n e^ n , n>2, 

and 

O O -1 

e = e 1 = 1. 

The algebra A cxt can be identified with C[r, r _1 ; e^, e^, . . . ]. 

Theorem 9.1. Under this identification, the operator T) ext : A cxt — > A cxt tafces i/je 
/orm 

where 
and 

with 
and 



d 2 d \ 1 

3ext = ( r ^ + ( c - r )7r +-®°» 

ar z or / r 



D° = 4° __ + — 

^ mn <9e° <9e° ^ n <9e r c 

m,n>2 m n n>l r ' 



min(m,n) — 1 

^mn = -rane^e° + (m + n-1- 2k)e 

k=0 



m+n—l—kk 



B° = -n(n - 1 + c)e° + (z - n + \){z' - n + l)e° n _ 



i- 



This result shows (at least on algebraic level) that the process X(t) is the skew 
product of the one-dimensional Laguerre diffusion with parameter c = zz' > and 
a Markov process on the Thoma simplex Q, which is generated by the operator T)°. 
This can be compared to the splitting of the multidimensional Brownian motion into 
the skew product of a one-dimensional diffusion (a Bessel process) and the spherical 
Brownian motion, see, e.g., [TM65J. 

In more detail: The algebra A ext is realized, in a natural way, as an algebra of 
functions on Q \ {0}. In this realization, elements of the subalgebra A° turn into 
homogeneous functions of degree with respect to homotheties of the cone Q. Since 
the Thoma simplex Q is a base of the Thoma cone, we may regard A° as an algebra 
of functions on Q. Note that A° is dense in the Banach space C(Q) of continuous 
functions on Q. Thus, 2)° becomes a densely defined operator in C(Q). As shown in 
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IB0Q9], the closure of £>° generates a diffusion process X°(t) in the Thoma simplex. 
□ Continuing the analogy with the multidimensional Brownian motion one can say 
that in this picture, X°(t) is a counterpart of the spherical Brownian motion and 
the variables e° play the role of spherical coordinates. 

10. Correlation functions 

Set R* = R\{0} and let Conf(R*) be the space of locally finite point configurations 
on R*. Following [BO0Q], we define a projection Q — > Conf(R*) by setting 

u = (a, 0, 5) h-> {on : on ^ 0} U {-ft : ft ^ 0}. 

Two points of Q are mapped to one and the same configuration if and only if they 
differ only by the value of the coordinate 5. Consequently, the restriction of the 
projection on the subset 

is injective. 

The equilibrium version X stat (t) of the process X(t) is obtained by taking the 
stationary distribution W as the initial one. The process X stat (t) is stationary 
in time; moreover, since W is a symmetrizing measure, one may extend the time 
parameter t from the half-line [0, +oo) to the whole real line R. We know that the 
stationary distribution W is concentrated on the subset Q'. It follows that the finite- 
dimensional distributions of the equilibrium process may be viewed as probability 
measures on the spaces (Q') k , k — 1,2, ... . Next, applying the above projection, 
we may interpret every k- dimensional distribution (k = 1,2,...) as a probability 
measure on 

(Conf(R*)) fc = Conf( R* U ■ ■ ■ UR* ), 

k times 

which is again a space of configurations. Such measures can be described in terms 
of the correlation functions. In other words, these are the dynamical or space-time 
correlation functions of a time-dependent point process. 

Theorem 10.1. The space-time correlation functions of the equilibrium process 
X sta,t (t) are determinantal. That is, they are given by minors of a kernel K(s, x; t, y), 
where s and t are time variables and ijeR*. 

The correlation kernel K(s,x;t,y) is called the extended Whittaker kernel. It 
was derived in |BO06at Theorem B] as a scaling limit of the correlation kernels 
of some equilibrium Markov jump processes on partitions. Explicit expressions for 
K (s, x; t, y) are contained in that paper. 

Note, however, that the paper [BO06aJ left open the question whether the kernel 
K(s, x; t, y) determines a Markov process (a subtlety here is that, in principle, it may 



In |BO09j . the operator 2)° is written down in a different coordinate system. 
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happen that the Markov property is destroyed in a limit transition). The results of 
the present section settle this question in the positive. 



11. The Meixner symmetric functions 

Let Z + C Z denote the set of nonnegative integers. Fix two parameters b and £, 
where b > as before and < £ < 1. 

The classical Meixner polynomials M n (x) are defined as the orthogonal polyno- 
mials corresponding to the following discrete probability measure supported by Z + : 

where 5 X denotes the Dirac measure at x (see, e.g., [KS96J). The measure w M is is 
known under the name of the negative binomial distribution. We use the standard- 
ization in which the M n 's are monic polynomials: M n (x) = x n + . . . . 
Consider the following second order difference operator on Z + : 

D M f(x) = ^f-f(x + 1) + JL-f{x - 1) 
£(b + x) + x 

- — c m 

(the factor 1 — £ in the denominator is introduced to simplify some formulas below). 
D M is formally symmetric with respect to the weight function and annihilates the 
constants. The Meixner polynomials are eigenfunctions of this operator, 

D M M n = -nM n . 

Moreover, they can be characterized as the only polynomial eigenfunctions of D M . 
The N-variate symmetric Meixner polynomials 

M v (xi, x N ) = M v \ NM (xi, . . . , x N ) 

are introduced following the recipe (13. ip . 
Set 

z +,ord = {0&i, • • • , xn) e : x\ > ■ ■ • > x N } 

and regard polynomials from as functions on Z^ ord . Then the Meixner polyno- 
mials M v become orthogonal polynomials with respect to the atomic measure 
on Z^ d defined according to (13.21) . 

The polynomials M v are eigenfunctions of an operator '■ A N — > A N , which is 
defined according to (13. 3p : 

D™M U = -\v\M v . 
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This operator can be realized as a difference operator on Z^ ord acting on a function 
/by 

N N 

D%f(x) = Mx)f(x + Si) + Bi{x)f{x - Si) - Ox) fix) 

i=l i=l 
N N 

= M*)(f(* + £ i) - f(*)) + E B i(*)(f(* - £ i) - /(*))■ 
j=i i=i 

Here x = (xi, . . . , x^v) ranges over Z^ ord , {ei, . . . , £#} is the canonical basis in M N , 
and the coefficients are given by 

£,(b + Xi) yj Xi—Xj + 1 



Mx) = n 



nn rf> nn \ 



1 J l-£ 2 

We will need a modified version of the truncation map ir^ : A — > Ajv; this is an 
algebra morphism n' N : A — >• Ajy-, which is defined on the generators pu (Newton 
power sums) by 

TV 

n'Mixt, ...,x N ) = ^[{xi -N + l) k - (-i + \) k \. 

i=i 

Since the right-hand side is a symmetric polynomial, the definition makes sense. It 
can be better understood in terms of the realization A C Fun(Y), see Remark 112.21 
below. 

Theorem 11.1 (cf. Theorem 14. 3[) . Let z, z' , and £ be complex parameters. For an 
arbitrary partition v = (1/1,1/2,...), there exists a unique element Wl u G A, which 
depends polynomially on z, z' and rationally on £, and such that for any natural 
number N > £(u) and any b > and £ G (0, 1), one has 



fal=N, z >=N +b -l) = M > 



v\N,b&- 

Definition 11.2. We call the elements Wl u the Meixner symmetric functions with 
parameters z, z', and £. 

Next, we need the Frobenius-Schur symmetric functions. These are some inho- 
mogeneous elements FS V G A indexed by arbitrary partitions v and such that 

FS U = S u + lower degree terms. 
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For their definition, properties, and explicit expressions, see |ORV03j . In particular, 
one disposes of a simple explicit expression of the Frobenius-Schur functions through 
the Schur functions. 

Theorem 11.3 (cf. Theorem 14.51) . The Meixner symmetric function DJl u with pa- 
rameters (z, z',£) are given by the following expansion in the Frobenius-Schur sym- 
metric functions: 

where 

C'(v,r,z,z>,Z) = (-l)H-H ( T 4 7 ) H " N m II (z + c(D))(z> + c(n))). 

Definition 11.4 (cf. Theorem I4.6p . The Meixner operator 2) M : A — > A with 
complex parameters (z, z',£) is defined in the basis {FS U } of the Frobenius-Schur 
functions by 

D M FS l/ = -\u\FSu + -j— 3"7 J2 ( z + C ( D ))(^' + <n))FS v \ n . 

* a&>- 

From this definition one sees that D M preserves the filtration of A and depends 
polynomially on (z, z') and rationally on £, with the only possible pole at £ = 1. The 
operator D M is uniquely determined by these properties together with the following 
one: If z = N = 1, 2, . . . , z' = N + b-1 with b > 0, and £ e (0, 1), then D M preserves 
the kernel of the map ii' N : A — > A N and the induced operator in A N coincides with 
the iV-variate Meixner operator with parameter b. 

Theorem 11.5 (cf. Theorem 14. 7p . The Meixner symmetric functions are eigenvec- 
tors of the Meixner operator 2) with the same values of parameters (z, z',£). More 
precisely, 

® M Wl u = -\u\Wl u . 

Definition 11.6 (cf. Definition 15.11) . For any fixed (z,z',^) E C 3 with £ ^ 1, 
introduce the formal moment functional ip M : A — > C by setting 

V> M (1) = 1, V(^ M /) = for any fe A. 

where the Meixner operator D M is taken with the same values of the parameters as 

Theorem 11.7 (cf. Theorem 15.21) . For any Young diagram v, 

r(FS u ) = (^-) H H(z + c(D))(z> + c(D)) ■ fe) 2 
^ nei> VII-/ 



LAGUERRE AND MEIXNER SYMMETRIC FUNCTIONS 



19 



This formula provides an alternative (equivalent) way of introducing the moment 
functional. 

Theorem 11.8 (cf. Theorem 15 .3f) . For any two Young diagrams v and \i, 

^ M (mujg = s vlt f l[(z + c(d))(z' + c(d)), 

where b V[li is Kronecker's delta. 

This result shows that the Meixner symmetric functions are pairwise orthogonal 
with respect to the inner product in the space A defined by 

(f,g) = (f,g) z , z >,C-=^ M (f9), f,geA. 

The inner product is nondegenerate provided that z and z' are not integers. 

Remark 11.9. The two characterizations of the Laguerre symmetric functions from 
Remark 15.41 extend . with obvious modifications, to the Meixner symmetric functions. 

12. The orthogonality measure for the Meixner symmetric 

functions 

To speak about the orthogonality measure we have first to find an appropriate 
realization of A as an algebra of functions on a space. In the context of the Laguerre 
symmetric functions that space was the Thoma cone Q. Now the relevant space is 
different: it is the countable set Y of Young diagrams. 

We will need the notion of modified Frobenius coordinates of a diagram A G Y. 
This is a double collection (a; b) = (ai, . . . , a d ; bi, . . . , b d ) of half-integers, where d 
stands for the number of diagonal boxes in A, a« = Aj — i + \ equals the number of 
boxes in the ith row of A plus one-half, and foj is the similar quantity for transposed 
diagram A'. For instance, if A = (3, 2, 2) then (a; b) = (2|, |; 2|, l|). 

Definition 12.1. Let A be the unital algebra of functions on Y generated by the 
functions of the form 

d 

P*(A) = Pk (a; -b) := - (-6,) fc ], k = 1, 2, ... , 

i=l 

where (a; b) is the collection of the modified Frobenius coordinates of a diagram 
A G Y. Elements of A are called polynomial functions on Y |KQ94] . 

Consider the generators Pi,P2, ■ ■ ■ of A (the Newton power sums). The assignment 
Pk l— >■ Pk{- ) extends by multiplicativity to an isomorphism A — > A and hence defines 
an embedding of A into the algebra Fun(Y) of functions on the set Y. This is the 
desired realization. 
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Remark 12.2. Now we can explain the origin of the map ir' N : A — > Ajy introduced 
in Section [TTJ Similarly to the realization A C Fun(Y), realize A^v as an algebra 
of functions on Y^v C Y, the subset of Young diagrams with at most N nonzero 
rows, by letting the arguments Xj of A^-variate symmetric polynomials to be equal 
to Xi + N — i, where A ranges over Y^, i = 1 . . . , N. Then tt' n is implemented by 
the natural map Fun(Y) — > Fun(YAr) assigning to a function on Y its restriction to 
Yiv- 

For a diagram A G Y, denote by A + the set of the boxes that can be appended to 
A. As before, A - is the set of the boxes that can be removed from A. By dim A we 
denote the number of standard tableaux of the shape A. 

Theorem 12.3. Under the realization A = A C Fun(Y), the Meixner operator 
D M : A — > A with parameters (z,z',£) is implemented by the following operator in 
Fun(Y), which will be denoted by the same symbol, 



= A ( x > °)[/( A u °) - /( A )] + E B ^ w( A \ °) - /(*)]• 



□ga+ dga- 

where 

Mx, □, = ^ + co)(,' + co) (l d 7 + ( yj A . □ , a + , 

* dga- 

C(A) = ^((l + fl|A| +£«>)■ 

Definition 12.4. Let z, z', and £ be complex parameters, |£| < 1. The associated 
complex measure on Y, called the (mixed) z-measure, is defined by 

m z ^{\) = (i - o zz ' U( z + c ( Q ))(^' + c (°)) ■ ^ |A| (ixir) 2 ' A e Y - 

□eA VII-/ 

One can prove that 

^M^(A) = 1. 

AeY 

These measures were introduced in |BU00j . and some closely related measures on 
the finite sets of Young diagrams with a fixed number of boxes appeared appeared 
earlier in [KOV93J. The measures M z z > £ are a special case of Okounkov's Schur 
measures [OkOlj . 

The next theorem relates the measures M ZjZ >£ to the formal moment functional 
ip M with the same parameters (see Definition lll.6p . 
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Theorem 12.5 (cf. Theorem 16 .3f) . Let (z,z') be admissible and < £ < 1. Then 
the mixed z-measure M Z)Z r % is a probability measure, all functions from A C Fun(Y) 
are integrable with respect to M ZrZ r^, and 

i> M U) = X>Wi^WA), V/ G A C Fun(Y). 

AeY 

Moreover, A is dense in the weight Hilbert space £ 2 (Y,M ZjZ r^). 

This implies that (under the above assumptions on the parameters) M 2 2 /^ serves 
as the orthogonality measure for the Meixner symmetric functions. 

Remark 12.6. The classical univariate Meixner polynomials are autodual in the 
sense that, in an appropriate standardization, they are symmetric with respect to 
transposition of the index and the argument, which both range over Z+. The similar 
autoduality property holds for the Meixner symmetric functions viewed as functions 
on Y under the realization A = A C Fun(Y). 

Indeed, under this realization, there is a simple expression for the functions FS^ ■ ) 
on Y: 

f |A|! dimA//i 
FS,(\)= I (|A|-H)! dimA ' 1 ^" , 
[0, otherwise, 

where A ranges over Y; see [ORV03] . Change the standardization of Wl u by setting 
where 

C>; z, z', := (^-) H ^ ]J(z + c(D))(z' + c(D)) 

is a normalizing factor. Then the above formula for FS^{\) combined with Theorem 
111.31 yields the following explicit expression for the function DJl' u (A): 





M 


I 






(M 




Ml)!' 






A*l)! 



dimis/fi dimA/// 



x 



n 



dim v dim A ^ (z + c(D))(z' + c(D)) 



Clearly, this expression is symmetric under v -h- A: 

wi' u (\) = m' x (v), p, AeY. 
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13. The Meixner jump process on the set of Young diagrams 

Return for a moment to the classical Meixner polynomials M n (x) and the asso- 
ciated difference operator D M on Z + . Let X± (t) denote the birth-death process 
X\A(t) on Z + whose jump rates are the coefficients of D M . That is, the rates of the 
jumps x — y x + 1 and x — > x — 1 are equal to (1 — + x) and (1 — 

respectively. This is a well-known instance of a birth-death process with linear jump 
rates. The negative binomial distribution w M is the stationary distribution of X^(t). 
The transition function of X± (t) can be expressed through the Meixner polynomials 
according to formula (12.11) . where one has to substitute 4> n = M n and w = w M . 

More generally, the coefficients A* and Bi of the operator (see Section [TT1) 
serve as the jump rates of a jump Markov process X^it) on the set Z+ ord . 

Even more generally, the following result holds (sec [BO06aJ). Assume (z, z') is 
admissible and < £ < 1. We know that then M ZjZ >^ is a probability measure. Its 
support suppM 2 2 /^ is the whole set Y if (z, z') belongs to the principal or comple- 
mentary series, or a proper subset of the form or {A : A' G Y^} if (z, z') belongs 
to the degenerate series. 

Theorem 13.1. Let (z, z') be admissible and < £ < 1. Then there exists a jump 
Markov process X^ z , ^(t) whose state space is supp M Z)Z >£ and whose jump rates are 
the coefficients A(X, □) and B(X, □) of the Meixner operator D M . The measure 
M ZyZ i£ is an invariant and symmetrizing measure for X^ z , ^(i) . 

The fact that the Meixner operator H M E is diagonalized in the basis Wl u of 
Meixner symmetric functions gives an expression for the transition function, which 
we state for the case of nondegenerate parameters. 

Theorem 13.2. Let (z, z') belongs to the principal or complementary series and 
< £ < 1. The transition function P(t; A, x) of the Markov process X[^ z , ^(t) on Y 
can be written in the form form 



Here A, x range over Y, the M v 's are viewed as functions on Y in accordance with 
the realization A C Fun(Y), and the explicit expression for tp M (VJt u VJt u ) is given in 
Theorem \11.8\ . 



As well known, the Meixner polynomials M n (x) are discrete analogs of the classical 
Laguerre polynomials L n (x). Namely, fix parameter b > and let £ f 1. Then one 
has the limit relation 




14. Approximation Meixner -> Laguerre 



lim(l - O n M n ((l - O x) = L n (x), 



(14.1) 
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where the scalar factor (1 — £) n is used to keep the coefficient of x n to be equal to 1. 

The limit relation (114. ip can be easily derived from the explicit expressions for 
the Laguerre and Meixner polynomials, see [KS96J. On the other hand, (114.11) can 
be explained by convergence of the weight functions: under the embedding 

Z+^R, x t-t(l-£)x, (14.2) 

the push-forward of the negative binomial distribution w M on Z + with parameters 
(&, £) converges, as £ 1 1, to the Gamma distribution (T(b))~ 1 x b ~ 1 e~ x dx on M. + . 

One more explanation can be given in terms of the univariate Meixner and La- 
guerre operators. In the same scaling limit regime (114. 2p . as the mesh of the lattice 
goes to 0, the Meixner difference operator turns into the Laguerre differential oper- 
ator. 

We are going to formulate similar statements in the infinite-dimensional context. 

Let G : A — » A be the operator multiplying every homogeneous element by its 
degree. In accordance with this, the operator (1 — £)~ G : A — > A, which appears in 
the next theorem, acts in the mth homogeneous component of A as multiplication 
by (1 — £) -m , for each m G Z + . 

The analog of ffTO) is 

Theorem 14.1. Let v be an arbitrary partition, DJl u G A be the corresponding 
Meixner symmetric function with parameters (z,z',£) G C 3 , £ 7^ 1, and £ u G A be 
the Laguerre symmetric function with parameters (z,z'). One has 

hm(l-0 M (l-0~ G ^ = ^ 

where convergence holds in the finite- dimensional subspace of A consisting of ele- 
ments of degree < \v\. 

This is a direct corollary of Theorems 14.51 and 111.31 

As a corollary of Theorem 114.11 combined with Theorems 14.71 and 111.51 one gets 
the following analog of the approximation of the univariate Laguerre differential 
operator by the Meixner difference operator: 

Theorem 14.2. Let D M : A — > A be the Meixner operator with parameters (z, z', £) G 
C 3 , £ ^ 1, and D : A — > A be the Laguerre operator with parameters (z,z'). One 
has 

lim(l - n- G o 2) M o (1 - n G = ©. 

Here we mean simple convergence on arbitrary elements / G A. Note that both 
the pre-limit and limit operators preserve the filtration of A, so that application of 
the both operators to a given / is contained in a fixed finite-dimensional subspace 
of A. 

For e > define the embedding 

t £ :Y^n, X = (a; 6) ^ {a, 0,5), 
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where (a; b) = (ai, . . . , a^; b±, . . . , bo) are the modified Frobenius coordinates of A G 
Y, by setting 



The image t e (Y) is a discrete subset of Q. As e \. 0, it becomes more and dense in 
O. This is the analog of the embedding (114.21) (e = 1 — £). 

The analog of the approximation of the Gamma distribution by the negative 
binomial distribution is 

Theorem 14.3 (cf. Section 5 in [BO00]). Let (z,z') be admissible and £ e (0,1). 
-As £ J 1, t/ie pushforward of the measure M ZjZ i£ under the embedding : Y — > 
weakly converges to the measure W z>z '- 

Recall that the weak topology on measures means convergence on bounded con- 
tinuous functions. Actually, more can be proved: convergence holds on any test 
function on Q, which is continuous and grows, as u = (a, (3, 5) — > oo not faster than 
a power of 5. In particular, as test functions one can take elements of A C Fun(u)); 
then, by virtue of Theorem I6.3l and ll2.51 the claim of the theorem means convergence 
of the formal moment functionals, 



which agrees with Theorems 15.31 and 111.81 

Finally, the approximation Meixner — > Laguerre holds on the level of Markov 
dynamics, which is used in the proof of the very existence of the diffusion process 
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